Dynamics of the collective modes of an inhomogeneous spin ensemble in a cavity 
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We study the excitation dynamics of an inhomogeneously broadened spin ensemble coupled to 
a single cavity mode. The collective excitations of the spin ensemble can be described in terms of 
generalized spin waves and, in the absence of the cavity, the free evolution of the spin ensemble 
can be described as a drift in the wave number without dispersion. In this article we show that 
the dynamics in the presence of coupling to the cavity mode can be described solely by a modified 
time evolution of the wave numbers. In particular, we show that collective excitations with a well- 
defined wave number pass without dispersion from negative to positive valued wave numbers without 
populating the zero wave number spin wave mode. The results are relevant for multi-mode collective 
quantum memories where qubits are encoded in different spin waves. 



PACS numbers: 37.30. +i, 42.50.Pq, 03.67.Lx 

I. INTRODUCTION 

The interaction between a single quantum system 
and a continuum of independent systems constitutes 
the standard model of dissipation and decoherence. 
The Weisskopf-Wigner model thus explains Markovian 
atomic decay due to coupling to the (oscillator) modes of 
the quantized radiation field, and the Caldeira-Leggett 
[1] and spin-star [2, 3] models investigate in detail the 
non-Markovian features due to the structure of the reser- 
voir density of states. Equivalent models have been used 
to describe the interaction of a two-state fermionic parti- 
cle with a bosonic continuum [4] , and the interaction of a 
laser mode with a gain medium [5-7] . Common to these 
works is that the focus is on the reduced dynamics of the 
central system, that is, on the decoherence and decay 
of the central system, the properties of the renormalized 
fermionic particle, or the behavior of the laser mode. 

While methods in quantum optics have been developed 
to address properties of the ensemble of radiation modes, 
such as the optical spectrum of fluorescence [8] , only few 
detailed studies of the ensemble dynamics exist. Recent 
proposals to manipulate and protect quantum informa- 
tion stored in collective degrees of freedom of ensembles 
urge a more precise description of the collective evolution. 
Systems of interest include the nuclear spins in a host 
material collectively coupled to a nitrogen-vacancy (NV) 
center [9] or a quantum dot [10], trapped atoms [11, 12] 
or polar molecules [13, 14] or electron spins [15, 16] cou- 
pled to a transmission line resonator, and ions in optical 
cavities [17]. 

The inhomogeneities of spins embedded in solids fol- 
low typically from properties of the host material, while 
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FIG. 1. (Color online) (a) The model studied in this paper 
describes a set of N+ 1 quantum oscillators with different fre- 
quencies interacting in a star topology. We assume that the 
width of the distribution of the outer oscillator frequencies 
can be controlled by some external parameter, (b) A physi- 
cal motivation for the model studied is a spin ensemble in a 
cavity exposed to a controllable linearly varying Zeeman shift 
provided by an adjustable magnetic field gradient [15]. 



very well isolated two-level systems may only have dif- 
ferent excitation energies if we employ a spatially de- 
pendent controllable perturbation to the system. Such a 
controllable inhomogeneity could be provided by, e.g., a 
magnetic field gradient. Spin systems with controllable 
inhomogeneities have been suggested as classical spin 
echo [18] and photon echo memories [19-21], and it has 
been proposed that the collective interaction with a cav- 
ity could be used to implement a holographic quantum 
register [14, 15]. In the latter, the qubits are stored as 
collective excitations equivalent to different spin waves, 
and a Raman process or a linear magnetic field gradient 
can be applied to deterministically change the wave num- 
ber between different collective spin wave modes. Such 
a memory can be read out by reversing the inhomogene- 
ity via applying an inverse Raman process or magnetic 
field gradient across the ensemble or, as pointed out and 
demonstrated in Ref. [22], by employing spin echo tech- 
niques to effectively invert the inhomogeneity. 

In this paper we consider the collective spin dynamics 
in an inhomogeneously broadened system of N indepen- 
dent spins Sj coupled to a single cavity mode as illus- 
trated in Fig. 1. In such a system a particular collective 
mode of the spin ensemble will experience an enhanced 
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coupling to the cavity mode due to constructive interfer- 
ence of the single spin coupling terms, and we will refer 
to this mode as the superradiant spin mode. The cou- 
pling to the cavity mode results in an energy gap between 
the superradiant mode and the other, subradiant, collec- 
tive modes of the spin system. A sufficiently large energy 
gap protects the superradiant mode from the dephasing 
caused by the inhomogeneous precession of the different 
spin ensemble members, an effect proposed in [10] for a 
nuclear spin ensemble interacting with a single electron 
spin. Consider an excitation created in the superradiant 
mode and left to evolve freely, i.e., without coupling to 
the cavity field, for a time T much larger than the inverse 
inhomogeneous width of the spin ensemble. If the cavity 
coupling is enabled at this time, the excitation will not 
interact with the cavity because of the destructive inter- 
ference of the individual spin components. Reversing the 
sign of the inhomogeneity or applying a 7r-pulse to the 
entire ensemble will, however, cause the system to refo- 
cus into the superradiant mode after another period of 
length T (spin echo). Note the significance of the cav- 
ity coupling during this process: If the cavity mode is 
not coupled to the spins, the excitation may refocus per- 
fectly into the superradiant mode where it would couple 
strongly to the cavity mode. If, however, the cavity cou- 
pling is turned on during the refocussing period, the en- 
ergy gap due to the coupling prevents the refocusing into 
the superradiant mode. Our main goal in the following 
is to understand the interplay between the spin-cavity 
gapping interaction and the inhomogeneity under such 
conditions. 

The paper is organized as follows. In Sec. II, we 
present the Hamiltonian describing our physical system 
and introduce the bare- and dressed-time states which 
arc related by an analytical expression derived in the Ap- 
pendix. In Sec. Ill we then consider the propagation of 
a spin wave packet and introduce a 'shear' approxima- 
tion to provide a concise quantitative explanation of the 
numerically observed effects. Section IV concludes the 
paper. 



II. PHYSICAL MODEL 

A large ensemble of spins prepared with high prob- 
ability in the spin down eigenstates is conveniently de- 
scribed by the Holstein-Primakoff approximation which 
replaces the spin operators by bosonic oscillator opera- 
tors, S^} w Oj, Sz^ — OjCij — 1/2, where a* is the cre- 
ation operator for a quantum oscillator representing the 
j-th spin [10, 23] and we set h = 1. This representation 
clearly does not reflect the saturation by a single exci- 
tation of each individual spin, but with a total number 
of excitations much smaller than the number of spins, 
the most significant state components will have at most 
a single excitation in each oscillator. In this limit, the 
Holstein-Primakoff approximation is a valid and conve- 
nient description of the spin ensemble as a collection of 



oscillators. 

We will assume that the spins are strongly coupled to 
a single mode of the cavity, and a cavity photon makes 
multiple round trips before it gets absorbed by the en- 
semble, so that we can neglect propagation effects in the 
medium. Wc further assume that the excitation energies 
of the spins are much larger than the strengths gj of the 
couplings to the cavity mode a c , so that in the rotating 
wave approximation the interaction picture Hamiltonian 
is W = H.Q+V where "H and V describe the free evolution 
and interaction, 

N 

Ha = uj c a\a c + ^^idjOjCLj, (la) 

N 

V = 5^ffjato c + H.c. =0(a c S t +H.c). (lb) 

Here the spin excitation energies ujj and the cavity mode 
frequency lo c arc defined relative to a common fixed fre- 
quency, and £1 = (X^ilffjl 2 ) 1 ^ 2 i s the collectively en- 
hanced effective coupling strength between the cavity and 
the superradiant mode described by the creation opera- 
tor [15] 

JV 

& = 2 a A ' where a i = n ' ^ 

The system Hamiltonian H is purely quadratic in the 
bosonic operators and does not describe any interac- 
tion between excitations. The dynamics of the non- 
interacting excitations is therefore completely described 
by the evolution of the one-excitation subspace, and to 
simplify the notation in what follows we will consider the 
evolution of states in that subspace. 

A. Bare-time states 

When the cavity is tuned far out of resonance and 
the spins precess freely in the presence of inhomogene- 
ity, the population in the superradiant state, defined as 
the singly excited state of the collective spin wave mode, 
|0) = 6*1 vac) gradually becomes subradiant. We will refer 
to the resulting distinguished subradiant states as bare- 
time states \t) and label them with the time r it takes 
for the system to reach them under this free evolution 
starting from the superradiant state, 

N 

\T)=U^(Tp\vac) = ^a je -^ r at|vac), (3) 
i=i 

where U^\t) = e~ in ° T is the free evolution operator. 
The bare-time states play an important role when the 
ensemble is used as a multi-mode quantum memory, and 
we note that if the spins are uniformly distributed along 
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the z axis and u>j — —KZj, the bare-time state |r) is the 
first excited state of the spin-wave with wave number 
k = kt [24]. In a finite size system, spin wave states, and 
hence the bare-time states, form an overcomplete basis in 
the set of singly excited spin states which are accessible 
from the superradiant state and 

N 

fain) =J2\ a tf e ~ iwATl ~ T2) = / e- MTl - r2 V ff H^, 

3=1 J 

(4) 

is proportional to the inverse Fourier transform of the 
effective spin density defined as the density of spin states 
weighted by the normalized coupling strength, 

N 

Poff(w) = ^2\aj\ 2 S(uj - ujj). (5) 

3=1 

In the following we shall be primarily interested in the 
limit of very many spins, and hence we shall indiscrimi- 
nately treat the spin excitation energies as an integration 
variable and a summation index. We will assume that the 
distribution of energies is such that the overlap (t2|ti), 
which describes how a collective excitation dephases due 
to the inhomogeneity, decays on a time scale T on the 
order of the inverse width of the ensemble. This in par- 
ticular implies that for |r| 3> T the bare-time state |r) is 
orthogonal to the superradiant state |0), the effect of the 
cavity on |r) is insignificant, and the time evolution of 
\t) is the same under the free Hamiltonian Hq as under 
the full Hamiltonian W + V. 

In the absence of coupling to the cavity, the time evo- 
lution of the bare-time states follows from the definition, 

U^(t)\r) = \t + r), (6) 

so the free evolution of the inhomogeneously broadened 
ensemble in the bare-time representation amounts to a 
displacement of the bare-time argument equivalent to a 
linear change of the collective spin wave number k <— > 
k + Kt. 

B. Dressed-time states 

To consider the ensemble dynamics under the full 
Hamiltonian, we introduce the dressed-time states |3>(t)) 
as the limit 

|$(r))= lim U(t)U^(-t)\r), (7) 

t— t oo 

where U(t) = e~ z ^ Ho+v ' t . In analogy to the time trans- 
lation property of the bare-time states as described by 
Eq. (6) we have for the dressed-time states 

U(t)Mr)) = m + T)), (8) 

so that in particular, once we have established the form 
of |$(0)), we can redefine the dressed-time states as 
|*(r))=W(r)|*(0)). 



si = o n = iott/t 




FIG. 2. (Color online) Time evolution of a collective spin 
excitation expanded in the bare-time (wave number) basis. 
The figures show the overlap |('r|^'(t))| 2 of the time evolv- 
ing spin state with the bare-time basis states. Figures (a) 
and (c) show that the evolution under the free Hamiltonian 
Ho (SI = 0) is a simple translation with time. Under the 
full Hamiltonian H (SI = 107r/T), the strongly coupled state 
jr = 0) is frozen, (b), while a bare-time state with initial con- 
dition | \E r (i — > — oo)) = \t = t) passes the r = region without 
ever populating the strongly coupled |r = 0) state, (d). Snap- 
shots of the expansion of the states in (c) and (d) are shown 
in Fig. 3. 



As the cavity field does not couple to bare-time states 
with large negative r, the bare- and dressed-time states 
coincide in the limit r — > — oo. In general, however, the 
dressed-time states are cavity polaritons with both spin 
and cavity components. As detailed in the Appendix, we 
can describe the coupling in terms of twice the complex 
phase of the inverse cavity mode susceptibility which can 
be computed as 

0H = 2arg^[ W - Wc -r2 2 ,5 c H] + i Sl 2 ^^j , (9) 

where we have introduced J c (oj) = 27rp G g(w) and 

<5cH = v\ P ^ J 1^', (10) 

J LO — U! 

with V denoting the principal part of the integral. In 
general <fi depends on the cavity frequency and coupling 
strength as well as the effective spin density, but for uj 
where ft 2 ^> \(ui — uj c )/(5 c (uj) + ij c (oj))\ we see that <f) is 
given by 

^( w ) = -2arctanQj5H) (11) 

independent of these parameters, and we will refer to 
the system as being in the 'strong coupling limit' if this 
condition is met for all w of interest. If the cavity is tuned 
far outside the support of p e g so that A c w lu c — u>, and 
|5 c (u;) + ij c (uj)\ is everywhere on the order of 2ir/T, we 
reach the strong coupling limit when SI 2 ^> \A c \2n/T. 
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In terms of 4>, we find in the Appendix that the Fourier 
transform of the dressed-time states can be expanded as 



l*M> = £ (e-^ l) |c) 



»dT + e- #(uj) / e luJT \T)dT. (12) 
Jo 



This expansion is readily evaluated for any choice of the 
spin couplings and density of states, and it can then be 
transformed back to provide the time dependent state of 
the system. Fig. 2 compares the time evolution of states 
in the absence and presence of the cavity coupling. The 
states are expanded on the bare-time (wave number) ba- 
sis, and we see that the free evolution is a simple transla- 
tion with respect to the bare-time argument. The cavity 
coupling, however, freezes the superradiant |t = 0) state 
and it significantly modifies the shape of a wave packet 
during propagation from negative to positive r, while 
asymptotically restoring its shape but causing an addi- 
tional translation of the bare-time (wave number) argu- 
ment. In Fig. 3, we show snapshots comparing the cou- 
pled and uncoupled evolution of wave packets starting at 
negative bare-time arguments, and we show that in addi- 
tion to the extra translational shear, the interaction pro- 
vides a change of sign of the wave packet. These numeri- 
cal results were obtained for an inhomogeneously broad- 
ened ensemble with p e g(uj) — — cos 2 (wT/2) l[_ xjr ](wT), 
see Fig. 4b. In the subsequent section, we will give an 
explanation of these observed phenomena and an analyt- 
ical approximation to the shear and phase shift in terms 
of the physical parameters of the model, and we will il- 
lustrate results for other inhomogeneity models. 



III. WAVE PACKET PROPAGATION 

To explain the results in Figs. 2 and 3, let us consider 
a spin wave packet with the asymptotic initial condition 



|*(f ->■ -co)) = / ip(t - t)|t) dr 



ip(t ;- t)|$(t)) dr. 

(13) 



When there is no coupling to the cavity field, e.g., be- 
cause the cavity is far detuned from the spin resonance, 
the spins evolve under the influence of only the free 
Hamiltonian, T-Lq, 



|*o(*)> - / i>{t-r)\r)dT, 



(14) 



as illustrated in Figs. 2a and c. Similarly, the evolution 
under the influence of the full Hamiltonian, H = Ho + V, 
can be expressed in terms of the dressed-time states as 



t = 4T 


/ At 




t = _ /\ 






t = -43/ 















r/T 



FIG. 3. (Color online) Evolution of a wave packet in the 
bare-time (wave number) basis in the presence and absence 
of a strongly coupled cavity. The initial state (bottom figure) 
is the negative bare-time state |r = — 4T) at t = — 4T, and 
the figures show the overlaps (r^^t)) (solid) and {r\^!o{t)) 
(dashed) at three different times during the evolution. With 
no cavity coupling, the spin wave packet evolution is a simple 
translation, while, the cavity coupling freezes the superradi- 
ant state, and unitarity forces the evolving wave packet to 
maintain a vanishing overlap (Oj >&(£)) = during the evolu- 
tion. The net effect on the final state (upper figure) is a phase 
shift 7r and an additional translation of the wave packet by 
At ~ T. 



When the cavity is strongly coupled to the ensem- 
ble, the cavity state |c) and the superradiant state |0) 
practically decouple from all other collective spin states. 
Therefore the time evolving wave packet (15) is prevented 
from refocusing into the superradiant state. We observe 
numerically (see Fig. 2) that the wave packet smoothly 
evades the superradiant mode and makes a transition 
from a superposition of negative to positive bare-time 
states. To a good approximation the wave packet asymp- 
totically regains its original shape up to a phase shift and 
a displacement of the expansion of the state on the bare- 
time basis states as illustrated in Fig. 3. 

To explain the observed shear effect and show how the 
displacement and the phase shift is related to the phys- 
ical parameters of the system, we will expand the spin 
component of the wave packet as J x(t, t)\r)dr. Using 
Eq. (12) we find the expansion coefficient 



X(r,t) 



= / V(*-r)|*(r))dr. 



(15) 



tp(t - t) for r < 0, 

i_ J e-iKt-r^+H^^p^duj for r > 0. 

(16) 

This is an exact description of the evolution of a wave 
packet in the spin ensemble, valid for an arbitrary time t 
and arbitrary wave packet envelope ip{t). In what fol- 
lows, we assume that ip(u>) has significant value only in a 
vicinity of uj within the support of p e g(u>), and that the 
linear approximation 

<f>(w) « 0(w o ) + Wo)0'(mo) = 0o + wAt (17) 
is valid in this vicinity. In this case, the integral over 
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frequencies in Eq. (16) can be approximately evaluated, 



B. Shear for delocalized wave packets 



X(r,t) 



%l>(t - t) 



for t < 0, 



e-^tpit + At - t) forr>0. 



(18) 



One readily sees that asymptotically, the wave packet is 
shifted by At and acquires a phase <fio, by the interaction 
with the cavity mode, and through Eqs. (9) and (17), 
these quantities are explicitly given by the physical pa- 
rameters of the model. 

In the following we will investigate the validity of 
the shear approximation (18) for three different exam- 
ple systems parametrized by their effective spin density 
Peg: (a) a uniformly broadened ensemble where the cou- 
pling strength of the individual spins is not correlated 
with the excitation energy of the spins, corresponding 
to p c ff(ui) = ^ ^(uoT) . (b) a uniform spatial dis- 
tribution of spins with constant magnetic field gradi- 
ent along a straight halfwave stripline resonator, corre- 
sponding to a trigonometric dependence on frequency, 
Peff(w) = ^ cos 2 (wT/2) l[_ 7rj7r ](u;T). (c) a spin ensem- 
ble with a Gaussian density of states, typical for a nu- 
clear spin ensemble in a semiconductor quantum well, 
where the inhomogeneous broadening is due to the hyper- 
fine interaction with a spin-polarized electron in the well 
(Knight shift). The frequency variance a 2 = (7r 2 -6)/3T 2 
is chosen to match the variance of p e s(u) for case b, the 
halfwave resonator. The properties of these three systems 
are summarized in Fig. 4. 



A. Shear for localized wave packets 

From the exact expression (16), we expect that the 
shear approximation (18) is valid if the frequency argu- 
ment of ip(cu) attains values within a narrow interval, cor- 
responding to states expanded on a wide range of bare- 
time states. To verify this expectation, we introduce the 
shear fidelity F as the asymptotic overlap of the exact 
asymptotic state \^(t)) with the bare-time states \t + At) 
maximized over the shift At, 



max lim I (t 

At t-HX> 



•Ai|*(t)>|. 



(19) 



As shown in Fig. 5, the size of the optimum shear indeed 
approaches <f)' as we consider wave packets which are well 
localized in frequency while at the same time the fidelity 
is seen to approach unity as the linear approximation to 
<j> becomes more accurate over the support of the wave 
function, exactly as we would expect. The configuration 
in which the shear approximation is studied in Fig. 5 
is well outside the strong coupling limit as most clearly 
seen from the fact that the form of <f> is very different 
from the strong coupling value (f2 = oo). This illustrates 
that high fidelity shear effects can be observed for any 
coupling strength. Also, the figure indicates that varying 
the cavity frequency and/or the coupling strength allows 
very different shear effects, even for a constant system 
geometry as described by p c g. 



We have verified above that the shear approxima- 
tion is indeed valid when ip(u>) only has support over 
a narrow frequency range where a linear approximation 
of the phase <f) is reasonable. The examples shown in 
Figs. 2 and 3, however, arc for states with well defined 
r, [tp(t) — S(t)] so that ip(to) is constant and extends to 
all values of co, invalidating any linear approximation of 
<f> as illustrated by Fig. 4. In this section our aim is to 
qualitatively investigate why for some systems we nev- 
ertheless observe a relatively high fidelity shear effect as 
suggested by Figs. 2 and 3 and illustrated quantitatively 
in Fig. 6. 

One way to understand the shear effect for delo- 
calized wave packets is to note that the expansion 
in Eq. (13) is not unique: Performing the convolu- 
tion in frequency space, we can write (r'\^o(t)) as 
/ e - 4a; (*— r )^(oj)p fiS (u))du}, making it clear that only the 
values of ip{uj) where p e s(uj) is nonzero are significant. 
Therefore, we have the freedom to consider only such ini- 
tial wave packets for which the support of "4>(oj) is within 
the support of p e g(cj). For example, if the support of 
p e ff(ui) is the interval [—w/T,w/T] (Figs. 4a and b), then 
instead of 4>(uj) — 1 corresponding to i[)(t) — S(t), we may 
think of the indicator function tp'(uj) = 1-[—k/t,w/t](u), 
whose inverse Fourier transform is ip'(t) = smc(irt /T)/T. 
This in particular implies that the linear approximation 
(17) need only be valid for values where p e ff(^) is nonzero, 
so that consulting Fig. 4 we would expect a high shear fi- 
delity for case b in the strong coupling limit, as is indeed 
observed in Fig. 6b. Not surprisingly, Fig. 6c shows that 
the shear fidelity for case c, where p e ff is qualitatively 
very similar to that of case b, is also high. 



IV. CONCLUSION 

To summarize, we have established a convenient de- 
scription of the evolution of an inhomogeneously broad- 
ened spin ensemble coupled to a central spin or a cavity 
mode. The dispersion in eigenfrequencies of the individ- 
ual spins leads to propagation, which is effectively de- 
scribed as a translation in a bare-time basis, equivalent 
to a spin wave number basis, and for a range of param- 
eters, the joint effect of the inhomogencity and the cou- 
pling of the cavity to a single superradiant mode can be 
understood in the same basis: The superradiant mode 
is gapped and hence becomes stationary (or a dressed 
state pair of the superradiant spin mode and the cavity 
mode form coupled, stationary eigenmodes of the system) 
even in the presence of inhomogeneities. Rather than di- 
agonalizing the problem and finding all eigenmodes of 
the system, we observe that the bare-time basis states 
with sufficiently large positive and negative time argu- 
ment evolve by a simple translation without dispersion in 
the bare-time basis argument. When the evolution causes 
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FIG. 4. (Color online) Examples of the phase <f> in the strong coupling limit for three different system geometries discussed in 
the text. For each part (a)-(c), the effective spin density p e ff is illustrated by the inset and the plot shows the value of 0oo(w) 
as given by Eq. (11) (solid lines), together with the linear approximation (17) corresponding to the shear observed for fully 
delocalized wave packets [ip(t) = S(t)] in the system as discussed in Sec. IIIB and Fig. 6. 
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FIG. 5. (Color online) The shear approximation for a local- 
ized wave packet, (a) The value of j c (filled) and S c (solid) 
as given by Eq. (10) for the spin system with uniform distri- 
bution of ui, p c g(u}) — j2ir {ujT), as also considered in 
Fig. 4a. (b) The phase 0(oj) given by Eq. (9) at uu c = and 
n = {27r/3,7r,oo}/T (dotted) and at uj c = -it /AT, f2 = tt/3T 
(solid). Dashed lines are linear approximations of <j>{tS) ac- 
cording to Eq. (17) in two different frequency domains around 
ojq, where the wave packets illustrated in part (c) of the fig- 
ure have significant value. The slope of the lines is given by 
the At that maximizes the fidelity (19). (c) The frequency 
domain wave packets tp(u), and the shear fidelity F for these 
wave packets. 



this argument to approach the value r = of the super- 
radiant mode, a transient nodal structure forms until the 
system again becomes a simple bare-time state, but with 
a changed phase and a shear, that we can compute from 
the physical parameters of the model. Since it is a lin- 
ear effect for highly polarized spin ensembles, the shear 
occurs not just for single excitations, as studied in this 



At/T= 0.499 
F= 0.9908 




FIG. 6. (Color online) Shear fidelity for delocalized wave 
packets. Plot shows the overlap of |r + t) with (i) the un- 
shifted state \t) (dashed), (ii) the state \t + At) correspond- 
ing to the shear approximation (dotted), and (iii) the actual 
asymptotic state compensated for 0o, e'" t ' \ 1 if(t)} (solid), for 
the same three spin systems (a)-(c) as described in Fig. 4 in 
the limit of t — > oo. For all three systems the wave packet 
traverses the r w region with little dispersion as discussed 
in Sec. IIIB and a time shift At in good agreement with the 
results of Sec. Ill as illustrated by Fig. 4. 



manuscript, but for any quantum or classical excitation 
of the initial state, provided that the ensemble remains 
strongly polarized. 

The result is relevant for the application in multi-mode 
quantum memories, as it allows an assessment of how pre- 
viously stored quantum states can be moved with high 
fidelity between register positions, while other quantum 
states are being read into or out of the superradiant 
modes. This is a very useful property for recent strategics 
for quantum computing with qubits encoded in different 
collective spin modes [14, 15]. Recent experiments with 
inhomogeneously broadened spin ensembles and stripline 
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cavities have verified the strong coupling to different col- 
lective spin systems [25, 26], and we imagine that the 
gapping and shear effects should be observable in these 
systems. 
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Appendix: Analytical description of the 
dressed-time states 



In the following, we will establish an analytical expres- 
sion for the expansion of the dressed-time states (7) on 
the cavity and bare-time states. 

We consider the evolution in terms of the resolvent 
operator, Q{z) = (z — H)^ 1 where z is a complex number. 
By considering the form of Q(z) close to the real axis we 
can compute the forward and backward propagators as 

g ± {w)= T i [ e lut Q{±t)U{t)dt= lim g{w±irj), 
J 77^0+ 

(A.l) 

where 0(i) is the Heaviside step function [27]. The prop- 
agators are in essence Fourier-Laplace transforms of the 
evolution operator U and are related to the Fourier trans- 
form of U by 

W(w) = J e lLOt U{t)dt = j [Q-(u) - G+(uj)] . (A.2) 

In terms of U we can write the dressed-time states, or 
rather their Fourier transform, 



l*H) = / < 



dte lut \$(t)) = lim e~ lUT U(uj)\-T). (A.3) 



To calculate a closed expression for |$(u>)) we partition 
the reachable part of the single-excitation subspace into 
the subspace spanned by |c) and that spanned by the 
bare-time states |t), and denote the orthogonal projec- 
tions onto these subspaces by C and S, respectively. 



1. Evolution of the cavity state 

We will first calculate the restriction of Q to the sub- 
space spanned by the cavity state |c). Since the sub- 
space is one dimensional, the restriction of any oper- 
ator is given by a complex number, CQ(z)C = CG c (z) 
where G c (z) = (c\Q(z)\c). As Q(z) by definition satisfies 
(z-H)G(z) = 1, we have that C = C{z-W){C+S)Q(z)C 
and = S{z - H)(C + S)Q(z)C. Noting that CH S = 0, 



it follows that 

G c (z) 



C 



z - CH C - CR{z)C 

= (z - uj c - Rci^yr 1 

where 1Z(z) is the level-shift operator, 

S 



n(z) = v + v- 



(A.4) 



(A.5) 



S{z-U)S 

and R c (z) = (c\lZ(z)\c). The level shift can be evaluated 
explicitly in the basis a^|vac), 



du), 



(A.6) 



and it is an analytic function of the complex variable z in 
the whole complex plane except for the support of p e s(u>) 
on the real axis, where R c (z) has a branch cut such that 
R c (z) does not tend to same value from below and from 
above the cut, 



R c± (lu) 



lim R c (u>±iri) 



Q 2 [J c (w)Ti7cM/2]. 



(A.7) 

where 5 C and -f c are given by Eq. (10). Nu- 
merically, R c ± is most easily calculated as 
-2mQ 2 J : t {@(±t)T- 1 {p eS }(t)}, where T denotes 
the Fourier transform. By Eq. (A.4) we can write 
G c ±(uj) explicitly as 



Gc±M = 



1 



(A.8) 



from which the exact form of (c|W(w)|c) follows from 
Eq. (A.2). We note that G c + is the susceptibility of the 
cavity mode, and that A c (w) = fl 2 5 c (uj) and r c (cj) = 
fi 2 7 c (u;) describe the frequency shift of the cavity and 
the decay from the cavity into the spin ensemble as dis- 
cussed in greater detail in Refs. [28-30]. 



2. Expansion of the dressed-time states 

By algebraic manipulations similar to those used to ar- 
rive at Eq. (A.4), we find that Gs(z) = SG(z)S is related 
to G c by 

g s (z) = Qf(z) + fi 2 G c (z)4 0) (z)|0)(0|4 0) (z), (A.9) 

where Q^\z) is the resolvent corresponding to the re- 
striction of % to the spin subspace, Qg*\z) = S/(z — 
SUS) [27]. Expressing the propagators in terms of the 

Greens functions as gfl(uj) = =R / e iu}t Q{±t)uf> \t)dt, 
we see that 



o|4°2m|-t 



=4/' 



,iu)(t— r) 



e(±t)(r\t)dt 



(A.10) 



in the limit of r — > oo, since we are assuming The matrix element is given by Eq. (A. 10), so that 
that the overlap (0|r) vanishes for large r. In 
other words, the backward propagator does not con- 
tribute to the interaction term. Since <S|$(w)) = 
lim^oo ie^' L ^ T [Qs+{^) — Gs-(w)]\—t), it follows from 
Eqs. (A. 9), and (A. 10), that 



(c|$(w)) = r!G c+ (w) 7c (w). (A.13) 



5|$H) = / e^ T [1 - iQ 2 e{T)G c+ {Lu)-f c {u)] \t)cLt. 

(A.11) 

To calculate the cavity component of the dressed-time 
states we note that CHqS = 0, so that it follows from 
the Dyson expansion Q(z) = £ (0) (z) + G{z)VG w (z) that 

Cg±(u)S = g c ±(u)Vg$l(u). Then 



Finally, to arrive at Eq. (12) we note that 



1 - i£l 2 G c+ (w)7 c (cj) 



G c -M' 



(A.14) 



e~ WT Cg ± (u)S\-r) 



\c)QG c± (uj)e- iUT (0\gfl(uj)\-r). (A.12) 



which, since G c _(cj) = G c +(w)*, we can write as e l 't'( u ) 
where </>(w) = — 2 arg[G c+ (w)]. 
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